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ABSTRACT 

The photonuclear production of vector mesons in ultraperipheral heavy ion collisions is 



CN , investigated within the QCD color dipole picture, with particular emphasis on the saturation 



model. The integrated cross section and the rapidity distribution for the AA — > V AA (V = 
^ I p,uj,(f), J/'^) process are computed and theoretical estimates for scattering on both light and 
^ ' heavy nuclei are given for energies of RHIC and LHC. A comparison with the recent STAR 
■ data on coherent production of p mesons is also presented. Furthermore, we calculate the 
photoproduction of vector mesons in proton-proton collisions at RHIC, Tevatron and LHC 
^ \ energies. 
O 

1 Introduction 

In ultraperipheral relativistic heavy-ion collisions (UPC's) the ions do not interact directly 
^ \ with each other and move essentially undisturbed along the beam direction. The only possible 
' interaction is due to the long range electromagnetic interaction and diffractive processes (For 
a review see, e. g. Ref. [1]). Due to the coherent action of all the protons in the nucleus, 
the electromagnetic field is very strong and the resulting flux of equivalent photons is large. A 
photon stemming from the electromagnetic field of one of the two colliding nuclei can penetrate 
into the other nucleus and interact with one or more of its hadrons, giving rise to photon-nucleus 
collisions to an energy region hitherto unexplored experimentally. For example, the interaction 
of quasi-real photons with protons has been studied extensively at the electron-proton collider at 
HERA, with ^/s = 300 GeV. The obtained jp center of mass energies extends up to W^p ~ 200 
GeV, an order of magnitude larger than those reached by fixed target experiments. Due to 
the larger number of photons coming from one of the colliding nuclei in heavy ion collisions, a 
similar and more detailed study will be possible in these collisions, with W^n reaching 950 GeV 
for the Large Hadron Collider (LHC) operating in its heavy ion mode. Similarly, the strong 
electromagnetic fields generated by high-energy protons allow us to study photon - nucleon 
processes in proton - proton interactions in a large kinematical range. These events can be 
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experimentally studied by selecting events with low multiplicity and small total transverse 
momentum. 

Over the past few years a comprehensive analysis of the heavy quark [2, 3, 4, 5, 6] and vector 
meson [3, 7, 8, 9, 10, 11] production in ultraperipheral heavy ion collisions was made considering 
different theoretical approaches. In particular, much effort has been devoted to obtain signa- 
tures of the QCD Pomeron in such processes [5, 6, 9], which can be used to constrain the QCD 
dynamics at high energies. Understanding the behavior of high energy hadron reactions from a 
fundamental perspective within QCD is an important goal of particle physics. In the late 1970s, 
Lipatov and collaborators [12] established the papers which form the core of our knowledge of 
Regge limit (high energy limit) of QCD. The physical effect that they describe is often referred 
to as the QCD Pomeron, with the evolution described by the BFKL equation. Since this equa- 
tion predicts that for s ^ oo the corresponding cross section rises with a power law of the 
energy, violating the Froissart bound, new dynamical effects associated with the unitarity cor- 
rections are expected to stop further growth of the cross sections [13]. This expectation can be 
easily understood: while for large transverse momentum /c^, the BFKL equation predicts that 
the mechanism g ^ gg populates the transverse space with a large number of small size gluons 
per unit of rapidity (the transverse size a gluon with momentum k± is proportional to l/k\), 
for small kj_ the produced gluons overlap and fusion processes, gg ^ g, are equally important. 
Considering this process, the rise of the gluon distribution with transverse momenta below a 
typical scale, energy dependent, called saturation scale Qsat is reduced, restoring the unitarity. 
It is important to emphasize that Golec-Biernat and Wusthoff [14] have shown that a saturation 
model is able to describe the DESY ep collider HERA data, in particular the transition from the 
perturbative to the nonperturbative photoproduction region (For improvements of this model 
see Refs. [15, 16]). Moreover, its parameter-free application to diffractive DIS has been also 
quite successful [14] as well as its extension to virtual Compton scattering [17], vector meson 
production [18, 19, 20], charm and longitudinal structure function [21, 22] and two-photon col- 
lisions [23]. However, once other approaches containing very distinct assumptions also describe 
the same set of HERA data, the description of the QCD Pomeron still is an open question 
and it deserves more detailed analyzes. An alternative way to constrain the QCD dynamics 
are studies on electron - nucleus interactions, since the saturation effects are amplified in a 
nuclear medium - oc A^^^. Several works have been made along these lines, studying the 
behavior of the observables in the kinematic region of the planned eA coUiders at HERA and 
RHIC (For recent reviews see, e.g., Refs. [24, 25]). In particular, recently we have analyzed 
the nuclear exclusive vector meson production considering two distinct theoretical scenarios 
[26]. Our results have demonstrated that the experimental analyzes of nuclear exclusive vector 
meson photoproduction in the future electron-nucleus colliders eRHIC and HERA-A could be 
useful to discriminate between the different theoretical scenarios, mainly if heavy nuclei are 
considered. It strongly motivates the extension of these studies for ultraperipheral heavy ion 
collisions which can be analyzed in the current and/or scheduled colliders. 

Recently, the STAR Collaboration released the first data on the cross section of the coherent 
p production in gold - gold ultraperipheral collisions at ^/s = 130 GeV [27], providing the first 
opportunity to check the basic features and main approximations of the distinct approaches 
describing nuclear vector meson photoproduction. The main aspect is that real photons have 
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a complicated nature. In a first approximation, the photon is a point-hke particle, although in 
field theory it may fiuctuate also into a fermion pair (See discussions in Rcfs. [1, 28]). In the case 
where there is a photon transition in a colorless antiquark- quark pair, the propagation of this 
colorless hadronic wave packet in a nuclear medium can be treated either in the hadronic basis 
as a result of Gribov's inelastic corrections or in QCD in terms of the partonic basis, which are 
complementary. Lets bricfiy discuss these two representations (For a detailed discussion sec Ref. 
[29]). The time scale characterizing the evolution of a qq wave packet can be estimated based on 
the uncertainty principle and Lorentz dilation. The lifetime of the photon fluctuation is given 
by tc — ^/{Q"^ + frty)-! where v is the photon energy, my is the mass of the fluctuation and 
is the photon virtuality. It is usually called coherence time. Using light-cone kinematics we can 
define the coherence length, which is given by l,- = tc- Moreover, one cannot decide whether a 
ground state V is produced or the next-excited state V, unless the process lasts longer than the 
inverse mass difference between these states. In the rest frame of the nucleus, this formation 
time is Lorentz dilated and is given hytf — 2u/ {my, —my) . Similarly, we can define a formation 
length given by If — tf. In the hadronic basis, the same process looks quite different. The 
incident photon may produce different states on a bound nucleon, the V meson ground state or 
an excited state. Those states propagate through the nucleus experiencing multiple-diagonal 
and off-diagonal diffractive interactions, and eventually the ground state is detected. According 
the quark-hadron duality, we expect that these two descriptions to be equivalent. However, 
as these two approaches have been used assuming different approximations, their comparison 
may provide a scale for the theoretical uncertainty involved. Furthermore, it is important 
to emphasize that at high photon energy u, both and If greatly exceed the nuclear radius 
Ra, which implies in the partonic basis that the transverse size of the qq pair do not change 
during the interaction with the target. This enables one to introduce the QCD dipole picture 
[30], where the process is factorized into the photon fluctuation in a qq pair and the dipole 
cross section. These aspects become the interpretation in the partonic basis more intuitive 
and straightforward than in the hadronic basis. Since the photoproduction of vector mesons 
in ultraperipheral heavy ion collisions have been analyzed in literature using only the hadronic 
basis, it motivates the study of this process using the partonic one. 

Other important point which motivates our analysis is that in the current studies of vector 
meson production in UPC's the extrapolations of the predictions for LHC energies in general as- 
sume a power-hke behavior for the 7p (A) — > Vp (A) cross sections. Despite the good agreement 
for the currently available energies, an extrapolation to higher energies of the experimental fits 
implies a large growth for the cross section which it would violate the unitarity at sufficiently 
high energies. Therefore, dynamical modifications associated to the unitarity corrections are 
also expected to be present in vector meson production [18, 19, 20]. As discussed before, these 
effects should be enhanced in nuclear processes [13, 31, 26]. 

In this paper the photonuclear production of vector mesons in UPC's is investigated within 
the QCD color dipole picture [30], with particular emphasis on the saturation model [14]. The 

exclusive vector meson production by real and virtual photons is an outstanding process provid- 
ing important information on the transition region from the soft dynamics (at low virtualities 
of the photon Q^) to the hard perturbative regime at high [32, 33] (For a recent review see 
Ref. [34]). In principle, a perturbative approach is only justified if a hard scale is present in the 
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process, e.g. the photon virtuahty and/or a large mass of the vector meson. For photoproduc- 
tion of hght mesons, such scale is not present and one has to rely on non-perturbative models. 
On some pQCD approaches, as the saturation model, even this soft process can be described, 
where the transition is set by the saturation scale. In the color dipole approach, the degrees 
of freedom are the photon (color dipole) and meson wavefunctions as well as the dipolc-nuclei 
cross section. Such an approach enables to include nuclear effects and parton saturation phe- 
nomenon. The latter one is characterized by the typical momentum (saturation) scale Qsat, 
which it has been constrained by experimental results in deep inelastic scattering (DIS) and 
diffractive DIS [14]. Here, we will use an extension of the phenomenological saturation model 
for nuclear targets [35]. This model reasonably describes the experimental data for the nuclear 
structure function and has been used to predict the nuclear inclusive and diffractive cross sec- 
tions for heavy quark photoproduction [36]. The nuclear saturation scale, QsA, provides the 
transition between the color transparency and the saturation regimes in the nuclear scattering. 
Concerning vector meson production, our starting point are the recent works in Refs. [20, 26], 
where this approach was applied for the proton and nuclear case. It is worth mentioning that 
although light meson photoproduction to be a soft process by definition, it is consistently de- 
scribed in the QCD color dipole picture whether there is a suitable model for the soft-hard 
transition, as occurring in the the saturation model. As a by product, we extend our analysis 
for the photoproduction of vector mesons in proton-proton collisions at RHIC, Tevatron and 
LHC energies. 

This paper is organized as follows. In the next section, we present a brief review of the 
ultraperipheral heavy ion collisions and the main formulae to describe the photon - hadron 

process in these reactions. Moreover, the photon spectrum in proton - proton collisions is 
discussed. In Section 3 we discuss the photoproduction of vector mesons in the QCD color 
dipole picture and the saturation model is shortly reviewed. In Section 4 we present our 
results for the integrated cross section and the rapidity distribution for the AA — > AAV and 
PP ~^ processes, where V — p,u!,(f), J/'if. Moreover, it is presented a comparison of our 
predictions with the STAR data for coherent p photoproduction on nuclei at energy ^/sNN = 130 
GeV and a discussion concerning related approaches used in its description is addressed. Finally, 
in the last section we present a summary of the main results and conclusions. 



2 Photonuclear vector meson production at UPC's 

In heavy ion collisions the large number of photons coming from one of the colliding nuclei will 
allow to study photoproduction, with energies W^n' reaching to 950 GeV for the LHC. The 
photonuclear cross sections are given by the convolution between the photon flux from one of 
the nuclei and the cross section for the scattering photon-nuclei. The final expression for the 
production of vector mesons in ultraperipheral heavy ion collisions is then given by, 

caa^aav (^\lSm?j = j duj — — — g^a-^va {^^a = Scu^^aSnn 
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where uo is the photon energy with oo^ 



m 



V 



/A'-firUp and a/Snn is the ion-ion c.m.s energy. 



The Lorentz factor for LHC is 7^, = 2930, giving the maximum c.m.s. energy W^a ~ 950 
GeV. In this process we have that the nuclei are not disrupted and the final state consists 
solely of the two nuclei and the vector meson decay products. Consequently, we have that the 
final state is characterized by a small number of centrally produced particles, with rapidity gaps 
separating the central final state from both beams. Moreover, due to the coherence requirement, 
the transverse momentum is limited to be smaller than p-r = \/2/Ra, where Ra is the nuclear 
radius. Therefore, these reactions can be studied experimentally by selecting events with low 
multiplicity and small total px- 

The photon flux is given by the Weizsacker- Williams method [1]. The flux from a charge Z 
nucleus a distance b away is 



[CO, 6^ 



Z^a, 



dui d% 



Kl in) + ^ Kl (r,) 



(2) 



where •jl is the Lorentz boost of a single beam and 77 = ub/jL, Kq^i^tj) are the modified Bessel 
functions. The requirement that photoproduction is not accompanied by hadronic interaction 
(ultraperipheral collision) can be done by restricting the impact parameter b to be larger than 
twice the nuclear radius, Ra — 1.2 A^^^ fm. Therefore, the total photon flux interacting with 
the target nucleus is given by Eq. (2) integrated over the transverse area of the target for all 
impact parameters subject to the constraint that the two nuclei do not interact hadronically. 
An analytic approximation for A A collisions can be obtained using as integration limit b > 2 Ra, 
producing 



dN {cu) 2 Z^ae 



diij 



r^Ko{f,)K,{n) + l- (Kl{n)-K, 



(3) 



where f] — 2uj RaIil- It is worth mentioning that the difference between the complete numeric 
and the analytical calculation presented above for the photon flux is less than 15 % for the 
most of the purposes [1]. 

In a similar way, vector meson production also occur when considering energetic protons in 
pp{p) colliders [11]. In this case the photon spectrum is given by [37], 
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(4) 



with the notation Q = 1 + [ (0.71 GeV')/Q^iJ and Q^^^ = c^VIt! {l-2u:I^JS^) ] ^ [i^liL?. 
The expression above is derived considering the Weizsacker- Williams method of virtual photons 
and using an elastic proton form factor (For more detail see Refs. [11, 37]). It is important 
to emphasize that the expression (4) is based on a heuristic approximation, which leads to an 
overestimation of the cross section at high energies ( 11% at = 1.3 TeV) in comparison 
with the more rigorous derivation of the photon spectrum for elastic scattering on protons 
derived in Ref. [38]. For a more detailed comparison among the different photon spectra see 
Ref. 
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3 Vector meson production in the color dipole approach 



Let us consider the scattering process 7p Vp in the QCD dipole approach, where V stands 
for both hght and heavy mesons. The scattering process can be seen in the target rest frame 
as a succession in time of three factorizable subprocesses: i) the photon fluctuates in a quark- 
antiquark pair (the dipole), ii) this color dipole interacts with the target and, iii) the pair 
converts into vector meson final state. Using as kinematic variables the 7*A'^ c.m.s. energy 
squared s = W^j^ = {p+qY, where p and q are the target and the photon momenta, respectively, 
the photon virtuality squared = —(f and the Bjorken variable x = Q'^/{W^^ + Q"^), the 
corresponding imaginary part of the amplitude at zero momentum transfer reads as [30] , 

Im A (7P ^Vp) = Y.ldz d'r ^l -,{z, r, Q') aX'^^ix, r) ^l^iz, r) , (5) 

h,h 

where ^^^(2;, r) and \E'^^(2;, r) are the light-cone wavefunctions of the photon and vector me- 
son, respectively. The quark and antiquark helicities are labeled by h and h and reference to 
the meson and photon helicities are implicitly understood. The variable r defines the relative 
transverse separation of the pair (dipole) and z {1 — z) is the longitudinal momentum frac- 
tions of the quark (antiquark). The basic blocks are the photon wavefunction, the meson 
wavefunction, ^j^l^ dipole-target cross section, cr^^^^. 

In the dipole formalism, the light-cone wavefunctions f^^^z, r) in the mixed representa- 
tion {z, r) are obtained through two dimensional Fourier transform of the momentum space 
light-cone wavefunctions, 'i/j^ ji{z, k) which can be completely determined using light cone per- 
turbation theory (see more details, e.g. in Refs. [19, 20, 32]). On the other hand, for vector 
mesons, the light-cone wavefunctions are not known in a systematic way and they are thus 
obtained through models. Here, we follows the analytically simple DGKP approach [40], which 
is found to describe in good agreement vector meson production as pointed out in Ref. [20]. 
In this particular approach, one assumes that the dependencies on r and z of the wavefunction 
are factorised, with a Gaussian dependence on r (For a detailed discussion sec Refs. [20, 26]). 

Finally, the imaginary part of the forward amplitude can be obtained by putting the expres- 
sions for photon and vector meson (DGKP) wavefunctions into Eq. (5). Moreover, summation 
over the quark/antiquark helicities and an average over the transverse polarization states of 
the photon should be taken into account. The transverse component (the longitudinal one does 
not contribute for photoproduction) is then written as [20, 40] 



ImAT{s, t = 0) = J J dz al(^ fv friz) exp 



—ujj. r 



X [^[z' + zf]K,{er) + ^K,{er)\aZ^\x,r) , (6) 
my my J 

with cr^ip^^* being the dipole-proton cross section in the nucleon case and the dipole-nucleus 
cross section for scattering on nuclei. In the photoproduction case, e — nif, where m/ is the 
quark mass of flavour /. The corresponding parameters for the vector mesons wavefunctions 
{mv, ujt-i fv-i etc) are presented in Table 1 of Ref. [26]. 
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In order to obtain the total cross section, we assume an exponential parameterization for 
the small \t\ behavior of the amplitude. After integration over \t\, the total cross section for 
vector meson production by real/virtual photons in the nucleon (proton) case reads as, 

, , [ImA{s, t = 0)]2 

where (3 is the ratio of real to imaginary part of the amplitude and By labels the slope parameter. 
The values considered for the slope parameter are taken from the parameterization used in Ref. 
[18]. For the p case, we have taken a different value in order to describe simultaneously HI and 
ZEUS photoproduction data. 

In addition, Eq. (6) represents only the leading imaginary part of the positive-signature 
amplitude, and its real part can be restored using dispersion relations TZeA — tan(7rA/2) Xm^. 
Thus, for the (3 parameter we have used the simple ansatz, 

with Aeflf = Acfr(W^7Ar, Q^) the effective power of the imaginary amplitude, which depends on 
both energy and photon virtuality. The correction coming from real part in photoproduction, 
where only transverse component contributes, is about 3% for light mesons and it reaches 
13% for J/^ at high energies. It is worth mentioning that a different computation of the /3 
parameter, as in Ref. [20], produces a larger effect even in the photoproduction case. An 
additional correction is still required for heavy mesons, like J/^f. Namely, skewedness effects 
which takes into account the off-forward features of the process (different transverse momenta 
of the exchanged gluons in the t-channel), are increasingly important in this case. Here, we 
follow the studies in Ref. [41], where the ratio of off- forward to forward gluon distributions 
reads as [41], 

22Aeff+3 r (Aeff + I) 
^^''^ ^ r (Aeff + 4) ' 

and we will multiply the total cross section by the factor i?^ for the heavy meson case. 

In the case of nuclear targets, By is dominated by the nuclear size, with B ~ i?^ and the 
non-forward differential cross section is dominated by the nuclear form factor, which is the 
Fourier transform of the nuclear density profile. Here we use the analytical approximation of 
the Woods-Saxon distribution as a hard sphere, with radius Ra-, convoluted with a Yukawa 
potential with range a — 0.7 fm. Thus, the nuclear form factor reads as [7], 



F{<1 = V^) = [sin(gi?^) - qRACOs{qRA) 

A q 



_l + aV. 

where po = 0.16 fm~^. Consequently, the photonuclear cross section is given by 

[ImAnnc{s,t = 0)Y 



(10) 



a(7A^l^>l)- ^^^^°7/''~"^^ (l + /3^) r dt\F{t)\\ (11) 
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with tmin — (^y/2<^)^, where cu is the photon energy. 



Having introduced the main expressions for computing vector meson production in the 
color dipole approach, in what follows we present a brief review of the saturation model and its 
extension for the scattering on nuclei targets. In the present work, we follow the quite successful 
saturation model [14], which interpolates between the small and large dipole configurations, 
providing color transparency behavior, a^ip ~ r^, as r ^ and constant behavior, aaip ~ Uo, at 
large dipole separations. The parameters of the model have been obtained from an adjustment 
to small X HERA data. The parameterization for the dipole cross section takes the eikonal-like 
form [14], 



a 



proton f ^ 
dip 



X, r ) = ao 



exp 



QL(5) = (|)' GeV^ (12) 



where the saturation scale Q^^t defines the onset of the saturation phenomenon, which de- 
pends on energy. The parameters, obtained from a fit to the small-x HERA data, are cxo = 
23.03 (29.12) mb, A = 0.288 (0.277) and xq = 3.04 ■ 10"^ (0.41 ■ 10"^) for a 3-flavor (4-flavor) 
analysis. An additional parameter is the effective light quark mass, rrif = 0.14 GeV, which 
plays the role of a regulator for the photoproduction {Q^ = 0) cross section. The charm quark 
mass is considered to be nic — 1-5 GeV. A smooth transition to the photoproduction limit is 
obtained via the scaling variable [14], x — [{Q^ + 4m|)/((5^ + W^jv)]- 

The saturation model is suitable in the region below x — 0.01 and the large x fimit needs 
still a consistent treatment. Making use of the dimensional-cutting rules, here we supplement 
the dipole cross section, Eq. (12), with a threshold factor (1 — x)"*''"% taking nthrcs = 5 for a 
3-flavor analysis and nthrcs = 7 for a 4-flavor one. This procedure ensures consistent description 
of heavy quark production at the flxcd target data [42] . 

Let us discuss the extension of the saturation model for the photon-nucleus interactions. 
Here, we follow the simple procedure proposed in Ref. [35], which consists of an extension to 
nuclei of the saturation model discussed above, using the Glauber-Gribov picture [43], without 
any new parameter. In this approach, the nuclear version is obtained replacing the dipole- 
nucleon cross section in Eq. (5) by the nuclear one. 



nucleus ^ ^ 
dip 



X, r'; A)^2jd'b[l- exp [-^T^(&) a^'-'ix, r')] } 



(13) 



where b is the impact parameter of the center of the dipole relative to the center of the nucleus 
and the integrand gives the total dipole-nucleus cross section for a fixed impact parameter. 
The nuclear profile function is labeled by T^ib), which will be obtained from a 3-parameter 
Fermi distribution for the nuclear density [44]. The above equation sums up all the multiple 
elastic rescattering diagrams of the qq pair and is justified for large coherence length, where 
the transverse separation r of partons in the multiparton Fock state of the photon becomes 
as good a conserved quantity as the angular momentum, i. e. the size of the pair r becomes 
eigenvalue of the scattering matrix. It is important to emphasize that for very small values of x, 
other diagrams beyond the multiple Pomeron exchange considered here should contribute {e.g. 
Pomeron loops) and a more general approach for the high density (saturation) regime must be 



8 



RHIC: AuAu — > AuAu + (V„ = p.co.c^.J/'P) rhic: SiSi— > SiSi + (V„ = p,(o,^,J/^) 




y 



Figure 1: The rapidity distribution for nuclear vector meson photoproduction on UPCs in AA 
reactions at RHIC energy (^/Snn = 0.2 TeV^. 

considered. However, we believe that this approach allows us to obtain lower limits of the high 
density effects. Therefore, the region of apphcability of this model should be at small values of x, 
i.e. large coherence length, and for not too high values of virtualities, where the implementation 
of the DGLAP evolution should be required. Consequently, the approach is quite suitable for 
the analysis of exclusive vector meson photoproduction in the kinematical range of the planned 
lepton-nucleus coUiders (eRHIC and HERA-A) as well as UPCs at RHIC and LHC colliders. 
It is noticeable that the energy dependence of the cross sections is strongly connected with the 
saturation scale Qsa{W^n). Namely, the saturation effects are larger whether the momentum 
scale is of order or larger than the correspondent size of the vector meson and the energy growth 
of the cross section is then slowed down. 



4 Results and discussions 

In this section we present the numerical calculation of the rapidity distribution and total cross 
section for the photoproduction of vector mesons in ultraperipheral heavy ion and proton - pro- 
ton collisions. Our main goal is to obtain estimates from the QCD saturation model for vector 
meson photoproduction in the kinematical range of the colliders RHIC, LHC and Tevatron. 
Furthermore, a discussion concerning the comparison of the present results with the currently 
available models [7, 8] based on vector meson dominance (VDM) and parton-hadron duality is 
also presented. 

Initially, let us consider the nuclear photoproduction of vector mesons in ultraperipheral 
heavy ion collisions at RHIC {^/s — 200 GeV) and LHC {^/s — 5500 GeV) energies for different 
nuclei. In Figs. 1 and 2 one presents the predictions for meson rapidity distributions and 
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LHC: PbPb — > PbPb + (V„ = p,co,^,J/'P) 



LHC: CaCa — > CaCa + (V„ = p,co,(|),J/'P) 




in Table 1 one shows the corresponding integrated cross section. For LHC the considered 
nuclei are lead (Pb) and calcium (Ca), whereas for RHIC one takes gold (Au) and sihcon (Si). 
The distributions for light and heavy vector mesons are placed in the same plot for sake of 
comparison among their orders of magnitude. The final state rapidity is determined by the 
simple relation y = ^ In , ^ = In — , that is the ratio between the photon energy uj and the 

longitudinal energy transfer in the laboratory frame. Consequently, the rapidity distribution 

reads as, 

da {AA AAV) da (AA AAV) 

where ^ is calculated using Eq. (1). Interchanging the photon emitter and target corresponds 
to a reflection around y = 0, with the total cross section being the sum of these two possibilities. 

Let us discuss the numerical results. The rapidity distribution is characterized by a plateau 
at mid-rapidity at both RHIC and LHC, even for the light mesons as the p production. Such 
plateau is more pronounced at LHC. This is in contrast with the theoretical predictions of Refs. 
[7, 8], where due to the matching between the photon spectrum and the photonuclear cross 
section at lower (RHIC) energies a clear double-peak structure appears. For LHC, that feature 
disappears since lower photon energies contribute in a smaller amount to the total distribution. 
This behavior is not present in our results mostly due to the inclusion of the threshold factor in 
the dipole-nucleus cross section, which enforces the correct behavior at energies near threshold 
region. This factor also prevents the sharp cut-off at the fragmentation (forward rapidities) 
region, which has appeared in recent calculation considering the J/^ production [11]. 

The integrated cross sections can be contrasted with the theoretical estimations using 
GVDM plus Glauber-Gribov approach of Refs. [8] as well as the estimation of Ref. [7], which 
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HEAVY ION 


J/^ (3097) 


0(1019) 


u (782) 


P (770) 


RHIC 


SiSi 


3.42 fib 


612 ijh 


764 nh 


6.74 mb 




AuAu 


476 ijh 


79 mb 


100 mb 


876 mb 


LHC 


CaCa 


436 ijh 


12 mb 


14 mb 


128 mb 




P1)P1) 


11.5 ml ) 


998 ml) 


1131 ml) 


10069 mb 



Table 1: The integrated cross section for nuclear vector mesons photoproduction at UPC's at 
RHIC and LHC energies. 

considers VDM plus a classical mechanical calculation for nuclear scattering and uses as input 
for the 7P ^ Vp reaction an extrapolation of the experimental DESY-HERA fits for meson 
photoproduction. Initially lets consider the latter approach (see Tab. Ill in Ref. [7]). At RHIC 
energy and Si nucleus, our results are about 20 % lower for p and cu, whereas gives a larger (p 
cross section and almost the same J/^ cross section. However, the situation changes for gold 
nucleus, where the present results are about 50 % larger than the estimates in Ref. [7]. At 
LHC energy and for Ca nucleus, our results gives higher cross section by a factor of order 10 
%, whereas for lead nucleus the factor reaches a factor 2 for light mesons and almost similar 
for the J/^ meson. Basically, the values are quite similar for hght nuclei. For heaviest nuclei, 
the results overestimate those ones in Ref. [7] when one considers hght mesons and become 
similar for the J/^ case. These features can be understood through the theoretical procedure 
when considering the nuclear scattering. 

Despite the QCD color dipolc model to be a different approach than that used in Ref. [7], 
we can qualitatively to understand the discrepancy on the results by looking at the 7 A — VA 
cross section in both models, which is the input for UPC's calculations in Eq. (1). For the 
dipole approach, this is given by Eq. (5) replacing p by the nucleus A. Qualitatively, one 
can roughly approximate the wavefunction as peaked at the dipole sizes around r ~ l/my — 
l/(2m/) in the photoproduction case. Namely, \|/^* oc 5 (r — l/mv)S{z — 1/2) and its 
normalization can be obtained from the electronic decay width constraint. This is supported 
by recent phenomenological studies on meson production within the dipole models [20, 26, 29]. 
Therefore, under this assumption the integration on longitudinal fraction z and dipole size r in 
Eq. (5) can be carried out. This procedure gives the following, up to a normalization factor. 



7 dipole 



t=0 oc 



a^''^-^s,r'^l/ml;A)]' 

(15) 



dt 167r 

for the nuclear photoproduction of vector mesons. The equivalent expression in the VDM 
approach is given by (See Eq. (13) in Ref. [7]), 

-^^1*-'^°^ 167r ' ^^'^ 
where the meson- nucleus total cross section has been labeled as atot {VA). In Ref. [7] the 
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following expression has been assumed for this cross section (See Eq. (12) in Ref. [7]) 

(Ttot {VA) = j d%{l- exp [-atot {Vp) Ta(6)] } , (17) 

which drives the behavior on energy and atomic number of the photonuclear cross section. As 
discussed in detail in Ref. [8] the above expression implies that in the black disk limit the 
total cross section becomes equal tiR\i which is the prediction of the classical mechanics. In 
contrast, a quantum mechanical approach implies that in that limit a tot — ^ttR^. In Ref. [8] 
the coherent p production in UPC was studied using the generalized vector dominance model 
and the quantum mechanical Gribov-Glauber approach, which implies that atot {VA) is given 
in the large coherence length limit by 

atot {VA) = 2 y rf^ft 1 1 - exp [-^a^* {Vp) T^(6)] | . (18) 

On the other hand, in the QCD color dipole picture we have 

^S?^"^ {x,rc^^)=2 Jd%[l- exp [-1 [x, r ^ ^) T^(6)] } , (19) 

where the dipole-proton cross section is given by Eq. (12) for the saturation model. In particu- 
lar, cr^[°*°'^ = (To — 23 — 26 mb in the case Qsat(^) ~ '^v'j while in the case where (5sat(^) ^ 
it presents the color transparency behavior a^ip — Q^g^^/irriy. In what follows, we present a 
theoretical and numerical comparison among the approaches. 

The dipole-cross section is somewhat similar to the meson-hadron cross section, since it gives 
the probability of scattering of a color dipole (a qq pair) off a hadron. Therefore, we can consider 
them at a level of similarity. Hence, the main difference between Eqs. (17) and (19) comes 
from the rescattering procedure. In the dipole approach, one has the standard Glauber-Gribov 
formahsm, whereas in the VDM approach used in Ref. [7] a classical mechanical version has 
been used. In the color transparency regime, characteristic at intermediate energies (or light 
nuclei) and/or heavy mesons, both the color dipole and the meson-hadron cross section are 
small. This implies the first terms in the rescatterings are the leading ones and Eqs. (17) 
and (19) give quite similar results. On the other hand, at very high energies (or heavy nuclei) 
and/or light mesons the black disk limit is reached. Thus, one obtains for the dipole approach 
a^^p — 2tiR\ and for the vector meson dominance model atot{Vp) = 7rR\. In this limit case 
the ratio cr^A.°^'^/^^A^ reaches a factor four. Similar arguments have been claimed in Ref. [8]. 

Now we compare our results with those ones in Ref. [8] , where the main focus is on the p 
and J/^ production. For p the predictions are computed only for RHIC energy y/s^N ~ 130 
GeV and we will consider it later on. We can anticipate that their results are closer to ours 
since a Glauber-Gribov approach is used in describing the scattering on nuclei. For J/\E' the 
theoretical approach for the photonuclear production was the coUinear QCD double logarithmic 
approximation, where the ^A — > J/^A cross section is directly proportional to the squared 
nuclear gluon density distribution [45]. There, it was considered an impulse approximation 
(no nuclear shadowing) and a leading twist shadowing version. The impulse approximation 
gives a larger cross section at central rapidity (about a factor 4 higher for Ca and a factor 6 
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for Pb), while at fragmentation region both approximations match each other at LHC energy 
for Ca and Pb nuclei. Our results arc closer to their impulse approximation, which suggests 
nuclear shadowing could be weak for J/\l' production. This feature can be easily tested in the 
first experimental measurements of coherent J/^ production on UPC's at LHC. Concerning 
the integrated cross section, they found 0.6 mb for Ca nucleus and 70 mb for Pb at LHC. Our 
results are 0.44 mb and 41.5 mb, respectively. Thus, our results are about 15 % lower for Ca 
and also 40 % lower for Pb. The difference between the predictions comes from mostly from the 
distinct QCD approaches considered used and the different photon flux in the UPC calculation. 

Recently, the STAR Collaboration published the first experimental measurement on the 
cross section of the coherent p production in gold - gold ultraperipheral collisions at \/i = 130 
GeV [27], providing the first opportunity to test the distinct approaches describing nuclear 
vector meson photoproduction. In what follows we compare our results with these experimental 
data and confront them with other theoretical predictions currently available. In Fig. 3- (a) 
we present the momentum transfer behavior for the photonuclear p production, which is the 
input for the UPC calculation. It is obtained by unfolding the integration over \t\ in Eq. (11). 
The dependence is proportional to the squared of the nuclear form factor in Eq. (10). The 
STAR Collaboration performed an exponential fit da^A-^pA/ dt oc e~^°"'=l*l for this reaction, 
with an approximate gold radius of Rau = V^B-^mc = 7.5 ± 2 fm, and obtained a forward cross 
section da^Ait = 0)/dt = 965 ± 140 ± 230 mb/GeV^. Accordingly, we have found the value 
daf^Ait = 0)/dt = 923 mb/GeV^, which is consistent with the experimental results. The usual 
dip occurs at \t\ ~ 0.015 GeV, which was also found in the calculation of Ref. [8]. The rapidity 
distribution of the coherent p production on UPC's is shown in Fig. 3-(b). The plateau at mid- 
rapidity remains, in contrast with the double-peaked structure appearing in the calculations of 
Ref. [8]. However, the values at central rapidities are similar and of order da{y ~ 0)/(iy ~ 100 
mb. 

The energy dependence of the cross section for the photoproduction of the meson p in gold- 
gold ultraperipheral collisions is presented in Fig 4. For comparison, the experimental data 
from STAR Collaboration at RHIC at ^/snn = 130 GeV [27] is also shown. Our theoretical 
estimations in the curves take into account the experimental cuts. The cut on the momentum 
transfer \t\ < 0.02 GeV^ slightly reduces the cross section by a few percents. In the range of 
rapidities \y\ < 3 (left plot), at energy ^/sNN = 130 GeV, we have found (Tsat(— 3 < y < 3) = 410 
mb in good agreement with the STAR measurement crsTAR(— 3 < ?/ < 3) = 370 ± 170 ± 80 
mb. For the cut \y\ < 1 (left plot), we have obtained crsat(— 1 < 2/ < 1) = 221 mb, whereas 
the STAR result is (7star(— 1 < y < 1) = 140 ± 60 ± 30 mb. In this case, our result is 
about 35 % higher than the central value of the STAR measurement. The values presented 
here are somewhat similar to the ones obtained in Ref [8], which uses the generalized vector 
dominance model (GVDM) and the Glauber- Gribov approach, including in addition the finite 
coherence length effects. As discussed before, our calculation on the photonuclear cross section 

— > VA considers large coherence length. This is not the case at low photon energies released 
at RHIC and the finite length effect could further suppress the cross section. We believe this 
suppression should not be strong, once they obtained (7gvdm(— 3 < y < — 490 mb and 
cgvdm(— 1 < y < 1) = 170 mb, respectively. 



13 




0.04 



(a) 



10^ 



< 10 



•10 



10" 



(b) 



Figure 3: (a) Momentum transfer dependence of the p meson photonuclear production at 
RHIC; (b) Rapidity distribution of coherent p meson production in gold-gold collisions at RHIC 
(\/snn = 130 GeV) for QCD dipole picture and saturation model. 



Finally, lets now consider the photoproduction of vector mesons in proton - proton collisions 
for RHIC (yi = 500 GeV), Tevatron (^s = I960 GeV) and LHC = 14000 GeV) energies 
(For a similar analysis of the photoproduction of heavy quarks in pp collisions see Ref. [46]). 
In Figs. 5 and 6, we present our predictions for the rapidity distributions and in the Table 2 
our results for the integrated cross section. The distribution on rapidity has similar features as 
for the nuclear case, presenting a plateau at mid-rapidity and without double-peaked structure. 
Concerning the integrated cross sections, for J/^> production a comparison with the results 
obtained in the Ref. [4] is possible. It is important to emphasize that in that reference a 
parameterization for the total 7p J /'^p cross section is used as input in the calculations. In 
this case we have that our predictions are similar for RHIC energies, being approximately 10% 
larger for Tevatron and LHC energies, which is expected since the saturation model describes 
reasonably well the HERA data. Very recently, predictions for the p and photoproduction 
in proton - proton collisions have been presented in Ref. [39]. We have that our results for 
production are 20 % smaller for RHIC energies but similar for LHC energy. In contrast, we have 
that our results for p production reach a factor 2 smaller than the predictions from [39] for all 
energies. This difference is due to the sizeable importance of the threshold factor. If this factor 
is not taken into account, the behavior of the cross section near threshold pthres = (^p + ^v) 
is overestimated. This gives additional contribution in an integrated cross section, mostly in 
the p case where the threshold scale stays at low energies yufhres ~ 1-6 GeV. Namely, for light 
mesons, the differential cross section da/dy receives a larger contribution from the low photon 
energies region which it is quite sensitive to the threshold corrections. 

Finally, lets present a short discuss on the experimental feasibility of the reactions consid- 
ered here. Although the vector meson photoproduction at AA or pp collisions to be a small 
fraction in comparison to the total hadronic cross sections, the experimental separation of these 
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Figure 4: Energy dependence of coherent p meson production in gold-gold in UPC's at RHIC 
(y/SNN = 130 GeV) in the QCD dipole picture and saturation model. Experimental data from 
STAR Collaboration [27]. 

reaction channels is possible. The rates are large even after taking into account the respective 
leptonic branching rations and/or acceptance estimates. The results presented in Tabs. 1 and 
2 show that the coherent photoproduction of light mesons are very high at RHIC and LHC 
For instance, the exclusive p production coming from these reactions reaches 10% of the total 
nucleus-nucleus cross section at RHIC, whereas corresponds to 50% of the PbPb total cross 
section at LHC. As photoproduction is an exclusive reaction, N + N N + N + V {N = p, A), 
the separation of the signal from hadronic background would be relatively clear. Namely, the 
characteristic features in photoproduction at UPC's are low pt meson spectra and a double 
rapidity gap pattern. Moreover, the detection (Roman pots) of the scattered nuclei (or protons) 
can be an additional useful feature. In hadroproduction, the spectra on transverse momentum 
of mesons are often peaked around meson mass, px ~ rnv- A experimental cut pt < 1 GeV 
would eliminate most part of the hadronic background. Hence, the rapidity cut would enter 
as an auxiliary separation mechanism. This procedure is specially powerful, since there will 
be rapidity gaps on both sides of the produced meson. For numerical estimates on these cuts 
procedures, we quote Refs. [11, 39]. 



5 Summary 



In summary, we have calculated the rapidity distribution and integrated cross sections of ex- 
clusive photonuclear production of vector mesons in ultraperipheral heavy ion collisions within 
the QCD color dipole picture, with particular emphasis on the saturation model. This is mo- 
tivated by the good agreement of this model in systematically describing the current data on 
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Figure 5: The rapidity distribution for vector meson photoproduction on pp or pp reactions at 
RHIC ('x/;^^ O.STeV; and Tevatron (^/S^ ^ 1.96 TeVj energies. 



Vm (mv)/COLLIDER 


J/* (3097) 


0(1019) 


ou (782) 


P (770) 


RHIC 


6.57 111) 


212.1 111) 


361.7 111) 


2.52 lib 


Tevatron 


24.64 nb 


476 nb 


646 nb 


4.84 fj,h 


LHC 


132 nb 


980 nb 


1.24 fjh 


9.75 /ib 



Table 2: The integrated cross section for the photoproduction of vector mesons in p + p{p) 
collisions at RHIC, Tevatron and LHC. 

vector meson photoproduction in scattering on protons and its reliable estimates for scattering 
on nuclei. The cross sections for the A + A — > A + A + V {V = p,u!,4>, J/'if) process were 
computed and theoretical estimates for scattering on both light and heavy nuclei are given 
for RHIC and LHC energies. The rates are very high, mostly for light mesons and at LHC 
energies. In particular, we compare our prediction for the coherent p meson production with 
RHIC data at y^Sjviv = 130 GeV. The corresponding results are in good agreement with the 
experimental results when considered the cuts on momentum transfer and on rapidity. We also 
have contrasted our results with the current models which consider vector dominance (VDM) 
or generalized vector meson dominance (GVDM) and Glauber-Grivov formalism for the nuclear 
scattering. We have pointed out the main sources of the discrepancies among those models and 
our estimates. In addition, the cross sections for the p + p{p) ^ p + p{p)V process are com- 
puted and theoretical estimates are given for the hadronic colliders RHIC and LHC in their pp 
mode and Tevatron. Finally, the experimental feasibility and signal separation on the reaction 
channels presented here are shortly discussed. Although the rates are lower than hadroproduc- 
tion, the coherent photoproduction signal would be clearly separated by applying a transverse 
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LHC: pp — > pp + (V„ = p,a),(^,J/^) 
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Figure 6: The rapidity distribution for vector meson photoproduction on pp reactions at LHC 
energy ^/SNN — 14TeV. 

momentum cut pr < 1 and two rapidity gaps in the final state. 
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